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We address the problem of barrier tunneling in the two-dimensional 75 lattice (dice lattice). In 
particular we focus on the low-energy, long-wavelength approximation for the Hamiltonian of the 
system, where the lattice can be described by a Dirac-like Hamiltonian associated with a pseudospin 
one. The enlarged pseudospin 5 = 1 (instead of S = 1/2 as for graphene) leads to an enhanced 
"super" Klein tunneling through rectangular electrostatic barriers. Our results are confirmed via 
numerical investigation of the tight-binding model of the lattice. For a uniform magnetic field, we 
discuss the Landau levels and we investigate the transparency of a rectangular magnetic barrier. We 
show that the latter can mainly be described by semiclassical orbits bending the particle trajectories, 
qualitatively similar as it is the case for graphene. This makes it possible to confine particles with 
magnetic barriers of sufficient width. 

PACS numbers: 73.23.Ad, 73.40.Gk, 73.23.-b, 73.22.Pr 



I. INTRODUCTION 

Nanostructures in two-dimensional electron systems 
are usually patterned by using gate electrodes creating 
adequate electrostatic (scalar) potential barriers. Typi- 
cal examples include quantum wires and quantum dots 
in semiconductor heterostructures. Yet this procedure is 
not applicable in the case of graphene, a single layer of 
carbon atoms arranged in a honeycomb lattice (HCL). 1,2 
At low energies, the HCL can be described by a rela- 
tivistic Dirac-Weyl Hamiltonian. 3 In this limit an electro- 
static potential barrier of arbitrary height and thickness 
is fully transparent for low-energy electrons at certain in- 
cident angles. 4 This effect — known as Klein tunneling 5 
- has recently been demonstrated experimentally us- 
ing graphene heterostructures. 6 On the other hand, suf- 
ficiently wide magnetic barriers yield zero transparency 
and therefore make it possible to confine electrons. ' 8 

In this article we investigate confining properties of 
the so-called dice- or 73-lattice 9 ' 10 , c -f- Fig- l a - This 
lattice is described by a Dirac-like Hamiltonian 11,12 - 
cf. (1) below — resembling the one describing the HCL 
and accordingly graphene, but with an enlarged pseu- 
dospin S — 1 11 compared to S — \ for graphene. 1 
It turns out that the presence of an integer pseudospin 
5 = 1 has striking consequences on the Klein tunneling: 
(i) Electrostatic barriers become even more transparent 
compared to the HCL case and (ii) there is a regime 
of super-Klein-tunncling with perfect transmission inde- 
pendent of the angle of incidence. The latter is related 



to a negative refraction index v = — 1 and allows for the 
realization of a perfect focusing system. 

The paper is organized as follows. First, in Sec. II, we 
introduce the low energy description of the 73-lattice and 
its implications on scattering at a barrier are presented. 
In Sec. Ill we solve the scattering problem of an electro- 
static barrier and discuss the effect of Klein tunneling. 
In Sec. IV we compare the analytical results with a nu- 
merical calculation using the recursive Green's function 
method. Section V is devoted to the general features of 
the 73-lattice in presence of a uniform magnetic field and 
Sec. VI deals with the scattering problem of a magnetic 
barrier. 



II. LOW ENERGY DESCRIPTION OF THE 
7^-LATTICE 

The 7~3 lattice [cf. Fig. la] contains three inequivalent 
sites per unit cell. Two of these lattice sites (generally 
called rim sites A and B) are threefold coordinated while 
the third site (referred to as hub H) is connected to six 
nearest neighbors. The energy spectrum of the 71s lat- 
tice consists of two electron-hole symmetric branches, in 
addition to a unique non-dispersive band at the charge 
neutrality point. 11 The latter roots in the lattice topol- 
ogy, which allows for insulating states with finite wave 
function amplitudes on the rim sites and vanishing am- 
plitudes on the hubs. The reciprocal lattice of T3 has 
a hexagonal first Brillouin zone [cf. Fig. lb] with the 
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FIG. 1: (a) The 75 -lattice in the a;-j/-plane, character- 
ized by translation vectors vi = (3/2; — v3/2)ao and V2 = 
(3/2; \/3/2)ao, with lattice constant ao. Rim sites A and B 
(open circles and squares, respectively) have a lower connec- 
tivity of 3 compared to 6 for the hub site (solid circles). The 
75 lattice can be viewed as two nested HCLs (emphasized by 
the two shaded hexagons), (b) The energy spectrum of the 
bulk 75 lattice. The first Brillouin zone is indicated by the 
dashed line on the flat band. 



electron-hole symmetric bands touching in the six cor- 
ners of the hexagon (Dirac points) . As in graphene, there 
are two inequivalent Dirac points, labeled K and K' . 

The low-energy Hamiltonian for the Ts lattice is ob- 
tained by a long-wavelength approximation to the tight- 
binding Hamiltonian. For a given K point it reads 11 



n = v F S- [-iV--A 



V 



(1) 



(in units where H = 1). Here v-p is the Fermi velocity, the 
operator V = (d x ,d y ) is the differential operator in the 
xy-plane, and we consider particles with negative charge 
q = — e < 0. Moreover, V(x, y) represents a scalar 
electrostatic potential (proportional to a unit matrix in 
spinor space), while A(x,y) is the vector potential re- 
lated to the z-component of a magnetic field. Pseudospin 
components 11 of S, 
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form a three-dimensional representation of SU(2) 
and satisfy angular momentum commutation relations 
[S m ,S n ] = ismniSi. However, contrary to Pauli matri- 
ces, they do not form a Clifford algebra, i.e., {S n , S m } 

2<5n )m l3. 

The 73-model (1) generalizes graphene by enlarging 
pseudospin. Each plane-wave momentum eigenstate |p) 
(at V = and A = 0) is pseudospin polarized with re- 
spect to the p-direction along which it experiences three 



possible eigenvalues S p = (+1,0,-1) . They translate 
into pseudo-Zeeman split energies E — v-p(+p, 0, — p) 
constituting three bands. The "non-magnetic" level 
S p = comprises the topologically protected dispersion- 
less band. 

In the following, we investigate the transmission prob- 
ability through electrostatic or magnetic rectangular bar- 
riers of constant height and thickness d which we assume 
as homogeneous along the y-direction as the archetype 
of a tunneling obstacle (c.f. Fig. 2). Accordingly, the 
three-component eigenspinors of (1) have the structure 



(3) 



which comprises the amplitudes on the three sublattices 
of the 73 lattice, -0A, "0H, and -0b- 

In order to obtain the transmission probability at non- 
zero V or A, we need to solve for the scattering problem 
by matching wave functions across the rectangular bar- 
rier, c.f. Fig. 2. Strictly sharp potential steps would cause 
large momentum scattering and so violate our assump- 
tion of close vicinity to one chosen JC-point. Therefore, 
as is customary, 4 we assume ultimately potential varia- 
tions as smooth on the length scale of the lattice constant 
ao but as sharp on the Fermi wavelength Af = 2ttvf/\E\ 
which is large at low energies \E\. 

To this end, let us derive the adequate matching con- 
ditions for S — 1 particles obeying (1) that turn out 
to be different from Schrodinger particles and different 
from S = 5 Dirac particles. In the latter case continuity 
of the wave function is required across a potential step 
(together with continuity of the spatial derivative in the 
Schrodinger case). These familiar conditions need to be 
altered when dealing with the Hamiltonian (1). As usual, 
we integrate the eigenvalue equation Hip = Eip over the 
small interval x £ [— e, e] along the x-direction and let e 
eventually go to zero. This yields 



■0H(e) 
■ V'B(e) 
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for non-diverging scalar or vector potentials, V and A. 
While ipn must be continuous, Eq. (4) only demands con- 
tinuity for ipA + -0b, allowing still for redistribution of 
occupation probability between A and B rim sites across 
potential steps. A similar integration of Hip = Eip along 
the y-direction yields continuity of ipn and of ipA — ips as 
matching conditions along the x-direction. 

A discontinuity in the wave function components V'a/b 
seems unexpected and an explanation in terms of physical 
quantities is helpful. Therefore, we calculate the proba- 
bility current 
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(see Appendix A), that satisfies a continuity equation, 
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FIG. 2: (a) Sketch of the scattering region with an elec- 
trostatic barrier of width d. (b) Wave vectors corresponding 
to the different plane- wave components in ipx, ipn, and V'ni; 
see Eqs. (9)— (11), for E > Vq. (c) Scattering trajectory in 
the special case E — Vo/2. Note that the opposite orienta- 
tion of q in the barrier region is due to propagation of holes 
[(E — Vb) < 0] in II, compared to electrons (E > 0) in I and 
III. 



for the probability density \ip\ 2 = \iPa\ 2 + IV'bI 2 + IV'hI 2 - 
From the expressions of j x and j y we can conclude, that 
the matching conditions correspond to a conservation of 
the probability current perpendicular to the barrier. On 
the other hand, parallel to the barrier the currents need 
not be equal on either side. 

More generally, for a barrier perpendicular to the vec- 
tor n = (cos(a), sin(a!)) T , with a € (— ? , ^), the current 



j ■ n = V2v F Re (^Ae" iQ + ^e ia )] 



(7) 



is conserved across the barrier. Therefore, only the lin- 
ear combination (-0Ae _1Q + ipB^ la ) is continuous across 
the barrier while {ipxe~ la — V'Be 1 "), in general, is discon- 
tinuous. 



III. ELECTROSTATIC BARRIER 

Let us consider the case of a rectangular electrostatic 
barrier 



V(x,y) = V e(x)@(d-x) 



(8) 



of width d and height Vq (cf. Fig. 2). Here Q(x) is the 
Heaviside step function. Due to momentum conservation 



parallel to the barrier, the wave function can be written 
as i/;(x)e lkyV [cf. Eq. (3)], with homogenous y-component 
k y of the wave vector. For ip( x ) we assume in- and out- 
going plane waves at x < 0, 



(9) 



with reflection amplitude r, that are solutions to (1) at 
energy E ^ 0. Consequently, the wave function under 
the barrier, < x < d , at non-zero V is of the form 
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while at x > d we consider out-going plane waves, 

V ^ I 



(10) 



(11) 



of transmission amplitude t. In (9)— (1 1) ip — 
arctan(fcj,/fc 2; ) is the direction of incoming and outgo- 
ing wave vectors. The barrier alters the wave vector 

q = (q x , ky) in its x-component, q x = y 1 (E - Vq) 2 - k 2 , 

and direction 8 = arctan(fcj / /g a .) . Further, we define 
s = sign(i?) and s' — sign(_E — Vq), and assume E ^ Vq; 
the special case E = Vq is treated separately below. 

By imposing the matching conditions (4) to ipi and 
at x — and to ipu and ipm at x — d we obtain a set 
of linear equations for the parameters r, t, a, b. Solving 
these equations yields 

r = -(l- 7 2 ) (l-e 2i ^) V- 1 

a = -27(l + ss'7)X> _1 

6 = 27(l-ss'7) e 2iqd V- 1 

t = -Ass'je-'^-^V- 1 (12) 

with 7 = cos^/cosS and T> — e 2lqd (1 — ss'7) 2 — 
(1 + ss'7) 2 . Finally, from Eqs. (12) we obtain the trans- 
mission probability 



T T3 = \t\ 2 = 



7 2 + 3 (l-7 2 ) 2 sin 2 ( g d) 



(13) 



This result coincides with the corresponding expression 
obtained recently for the line-centered square (Lieb) lat- 
tice. 13 The functional form (13) differs considerably from 
the graphene case for the same rectangular barrier, 1 
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FIG. 3: Polar plot of the transmission probability T as a 
function of the injection angle tp for different values of the 
energy E for (a) the T3 lattice and for (b) the HCL lattice. 
The width of the barrier is dVo — 60 in both panels. The 
red dashed line is for E = O.lVb, the blue dotted line is for 
E = 0.4Vo, the green solid line is for E — 0.5Vb, and the 
black dot-dashed line is for E = 0.95Vo- For the Tz lattice at 
energy E — 0.5Vb we get perfect transmission, independent 
of the incident angle tp. 



The transmission Tj- 3 is plotted in Fig. 3(a) versus the 
incident angle tp for various energies < E < Vq , i.e., 
in the regime of ss' = — 1. Figure 3(b) shows the result 
(14) for the HCL for the same parameters for compari- 
son. For the 73-case, we find almost perfect transmission 
for a wide range of incident angles. For perpendicular 
incidence, i.e., tp = which implies 9 — 0, we always ob- 
serve perfect transmission T = 1 that is usually referred 
to as 'Klein tunneling'. The condition T = 1 is also ob- 
tained for resonant values of qd being equal to an integer 
multiple of n, as in graphene. Comparison between the 
two cases shows that electrostatic barriers on T3 are sig- 
nificantly more transparent than on graphene. We can 
even establish the inequality 



T T 3 > T, 



HCL 



for E G 



2 



holding at all tp for same barrier height Vq and thickness 
d. We conjecture that these high transparencies make 
T3 even "more ballistic" than graphene and less sensitive 
to disorder 14 in its transport properties at zero magnetic 
field. 

The most striking feature in the transmission proba- 
bility is that at energy E — Vq/2, where 7 = 1, we obtain 
perfect transmission T = I for any angle tp of incidence 
and any barrier thickness. This physical effect can be 
analyzed by considering the scattering of a plane wave at 
a single electrostatic interface V(x,y) = VoO(x) which 
describes a Tip-junction. At E = Vq/2 electrons propa- 



gating at x < and incoming angle tp continue to prop- 
agate as holes at x > and a diffraction angle 9 = tt — tp 
[cf. Fig. 2(c)]. Therefore, the system can be characterized 
by a negative diffraction index n: 15 



s s 



, sin tp 
sin 6* 



I . 



(15) 



This holds true for both, the T3 and the HCL lattice. 
However, for the HCL the reflection probability keeps an 
angular dependence 

Rucl = ^[1 - cos(2<p)} 

while in the case of the T3 lattice Rj- 3 = 0. 

Finally, we discuss the special case of E = Vq. Here let 
us first consider k y = 0, that is, tp = 0. Then the wave 
function in the barrier region is constant, 

*ii(x) = {a,h,b) T . 



(16) 



The requirement of current conservation yields the four 
equations 



2(1 + r) = a 
v^s(l-r) 



-b, 
h, 



a - 
h 



b = 2te 
V2 ste 



i/>- , d 



ik x d 



(17) 
(18) 



and 
= 1. 



which readily give r = 0, (a + b) = 2, h = \/2s 
t = c~ lfejjd and therefore perfect transmission, \t\ 2 
Note that (a — b), describing a finite current in the y- 
direction in the barrier region, a priori cannot be deter- 
mined. We may even include an arbitrary linear combi- 
nation of topological states, 



*n P (*) = £« 





V-v 



(19) 



which automatically satisfies ?/>a + V'b = for any q and 
any amplitudes a q . Due to the zero hub component, 
none of these states or any linear combination thereof 
carry current and therefore (19) does not participate in 
the scattering problem. 

At E = Vq and finite momentum k y , that is, tp =/= 0, 
the wave function ^u(x)e lkyV in the barrier region is a 
superposition of two evanescent modes and the topolog- 
ical states, given by 
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(20) 



where 9 q = arctan(fc y /g). As argued above, the vanishing 
hub component directly implies that this state does not 
carry current, hence we have perfect reflection (r = — 1) 
and zero transmission t — 0. Note that continuity of 
(tpA + iPb) at x = and x = d only makes it possible 
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to determine two of the (in principle) infinite amount of 
parameters a, 6, {a g }. We stress the particularity that 
at E = Vq, the transmission T(ip) only is non-zero (and 
equal to 1) at the singular angle cf> = 0. At E = Vq 
there is no tunneling contribution to the current at finite 
incident angles, which is in striking contrast to the case 
of graphene. 1 

Finally, we confirm to find the same transmission prob- 
ability (13) when considering a potential step along an 
arbitrary direction, perpendicular to some vector n, as 
expected from the rotational invariance of the low-energy 
Hamiltonian (1). 



IV. ELECTROSTATIC BARRIER IN THE 
TIGHT-BINDING MODEL 

So far we have assumed that the electrostatic poten- 
tial does not induce scattering between the two valleys, 
that is, between the two inequivalent K and K' points 
of the first Brillouin zone. In order to examine the va- 
lidity of this approximation and to study the effect of 
a smooth variation of the potential on the transmission, 
we also compute numerically the transmission probabil- 
ity through an electrostatic barrier in the tight-binding 
model. 

In the following we consider the effect of a potential 
V(n ■ r) on the tight-binding lattice, where n is an unit 
vector in the direction where the potential changes. Al- 
though such a potential is in principle translationally in- 
variant along the direction perpendicular to n, the dis- 
creteness of the lattice reduces the translational invari- 
ance to a discrete lattice symmetry. In particular, if n 
is commensurate with the discrete lattice, the combined 
system of potential and lattice is translationally invari- 
ant under shifts of t, where t is lattice vector perpen- 
dicular to n. The wave functions are then of the Bloch 
form ijj{r) = e 1 ' £ ii t ' , "/l t lu(r) with u(r + i) = u(r) and 
G [ — ]fr> iff )■ The original problem of an infinitely ex- 
tended system can thus be reduced to a finite unit cell of 
width |i| with periodic boundary conditions in the direc- 
tion of t. The problem must then be solved for each fc|| 
individually, with hoppings across the unit cell bound- 
aries being multiplied by a phase factor of e ±lfe nl*l. 16 As 
long as the Fermi momentum kp in the leads (measured 
with respect to the K and if' -points) is smaller than j|r , 
the scattering states have conserved parallel momentum 
in the leads that allows us to compute the transmission 
probability T as a function of the incident angle ip as in 
the continuum case. The transmission itself is computed 
using the numerical method of Rcf. [17]. 

We consider two orientations of the potential step: A 
step along the y-direction [t = Vx — v 2 with v 12 the 
lattice vectors from Fig. 1(a)] and along the x-direction 
(t = 2vi + v 2 ) which we — borrowing from the HCL 
nomenclature — call "zigzag" and "armchair" directions, 
respectively. Those two geometries, together with the ef- 
fective unit cells, are depicted in Figs. 4(a) and 4(b). 



(a) (b) 




FIG. 4: Comparison of numerical transport simulations and 
the analytical result for the continuum limit. Panels (a) and 
(b) show the geometries used in the numerical simulations, 
with the potential step along either (a) a zigzag direction or 
(b) an armchair direction. Panels (c)-(f ) show a comparison 
of the transmission T as a function of the incident angle tp 
as computed numerically from the tight-binding model (black 
squares, armchair; blue crosses, zigzag) and from the contin- 
uum result (13) [red line, left out for clarity in (f), where 
there is only perfect forward scattering at tp = 0]. The pa- 
rameters used are E = O.Olt, d = 1200^0, d s = VSa , 
and (c) V = 0.02i, (d) Vo = 0.018t, (e) Vo = 0.015t, and (f) 
V = O.Olt. 



While we only show results for two particular directions, 
we have checked that orientations in between give equiv- 
alent results. 

The two orientations zigzag and armchair are funda- 
mentally different on the lattice level: A step along a 
zigzag direction cannot scatter the valleys, as the K and 
if'-points are at different transverse (Bloch) momenta 
fen. In contrast, for a step along an armchair direction 
both K and if'-point are projected to fc|| = 0, hence in- 
tcrvalley scattering is always possible. Note that while 
an abrupt step in the zigzag direction cannot scatter the 
valleys, there still may be corrections to the low-energy 
Hamiltonian (1) in the form of mass-like terms, as an 
abrupt step makes the lattice sites inside a unit cell in- 
equivalent (for an extensive discussion of equivalent ef- 
fects in the HCL, see Ref. [18]). 

In order to compare our result (13) from the low-energy 
Hamiltonian to the lattice results, we restrict ourselves to 
small energies E <^ t and barrier heights Vq <C t, where 
t is the hopping matrix element of the tight-binding lat- 
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FIG. 5: (a) Transmission as a function of incident angle ip 
through a smoothed potential step along a zigzag (black line) 
and armchair direction (red line). Parameters are E = O.lt, 
Vo = 0.2t, d = 200\/3ao, and d B = 5.5V3a . (b) Probability 
for a valley flip in the scattering process (given by the to- 
tal probability of transmission or reflection into the opposite 
valley, Ti ntC r + Winter) for a smoothed potential step along the 
armchair direction at the incident angle if — 70° as a function 
of E. The remaining parameters are as in (a). 



tice. In addition, instead of abrupt jumps we ramp the 
potential linearly over a distance d s with kpd s < 1. A 
comparison between our analytical and numerical results 
is shown in Figs. 4(c)-4(f). We find perfect agreement 
between the low-energy, continuum and lattice results. 
As predicted from the rotational invariance of the low- 
energy Hamiltonian, the transmission probability is in- 
dependent from the direction of the potential step. In 
particular, we find a virtually angle-independent trans- 
mission for E = 0.5Vb; small deviations are only visible 
for angles close to grazing incidents, where lattice effects 
become important. We also find strictly forward trans- 
mission in the case of E = Vq as predicted by the low- 
energy theory. 

Barrier transmission on the T3 -lattice differs not only 
at low energies from graphene, as seen from the different 
functional forms of the transmission probabilities T(ip) 
[Eqs. (13) and (14)], but pronounced dissimilarities arise 
also at higher energies as investigated exemplarily now. 

While wc found that the transmission probability was 
independent of the lattice direction in the case of a small 
potential step at low energies, we find a quite pronounced 
angle dependence at higher energies and larger poten- 
tial steps or larger d s , as seen in Fig. 5a for parameters 
E = Vq/2 where the continuum theory yields unit trans- 
mission probability. Furthermore, for a potential step in 
the armchair direction the transmission is considerably 
smaller than for a step along the zigzag direction. 

We can ascribe the difference to intervalley-scattering: 
A step in the zigzag direction still excludes intervalley 
scattering due to translational symmetry but it may take 
place at a step in armchair direction. 

In Fig. 5b we investigate intervalley scattering in more 
detail for an armchair barrier and show the probability 
for a valley flip in the scattering process, that is, the prob- 
ability of transmission or reflection into the other valley. 
For energies E in the range < E < Vo pronounced res- 
onances are seen where intervalley scattering is strongly 
enhanced. This effect even increases with increasing d s 



(not shown), in stark contrast to the case of a HCL lat- 
tice, where smoother potential barriers reduce interval- 
ley scattering. We attribute this effect as being mediated 
by crossing of the Fermi energy with the flat band con- 
necting both valleys when V{r) = E for some r and 
this mediation taking place more efficiently at smoother 
barrier slopes (similar as in a Landau-Zener transition). 
For energies E > Vq, when this crossing does not take 
place anymore, intervalley scattering is suppressed, ac- 
cordingly. 



V. UNIFORM MAGNETIC FIELD 

Now we turn to the effects of a perpendicular mag- 
netic field B z (x,y) on 73. In the homogeneous case, 
B z (x,y) — Bo, Landau levels (LLs) appear, 11 which re- 
semble the relativistic Landau splitting of graphene in 
their proportionality to uj c cx VBq, where uj c = ^/2vf/(-b 
is the cyclotron frequency and £b = y/ c/(eBo) the mag- 
netic length. However, the non-zero Landau energies for 
the T3 lattice are 



ET 3 = ±UJr 



- for n > 1 
2 



(21) 



and differ crucially from their values of graphene given 
by E^ CL — ±uj c y/n ■ We mention that this latter spec- 
trum is basically known from massless Dirac-Fermions 19 
moving at the velocity of light; it can be attributed to 
a Zeeman-contribution ±ec.Bo for the two spin helicities, 
according to which the relativistic LLs at positive ener- 



2 i 2 ■ As a re- 



gies take the values E{ » j 
suit, each LL exhibits a doubled Zeeman-degeneracy with 
the exception of the zero-energy level Eo which remains 
non-degenerate. 20 The similar mathematical structure of 
the Dirac-Weyl equation for graphene causes E n — E^ ch 
together with only one pseudospin projection in the ze- 
roth LL which is responsible for the anomalous Hall con- 
ductivity 20 of graphene, as observed in experiment. 21 On 
the other hand, the spectrum (21) arises due to the ab- 
sence of a pseudo-Zeeman contribution in the T3 case. It 
is a signature of the pseudospin 1 physics 22-25 and has 
not yet appeared in other physical systems. 

In Landau gauge, A = (—B y, 0, 0) T , the LLs of T3 are 
given as 



( y/i. 



V2(2n-1) 



1 



</>n = sign(£:„)V2n - 1 4> n -i , n > 1 

n 4> n 



J 



where the <p n are eigenstates of the harmonic oscillator 
at frequency uj c and we define tf>—i = 0. While one half 
of the probability density resides on the hub sites, as ex- 
pected from the lattice structure, the remaining half is 
shared unequally among the lattice sites A and B accord- 
ing to (22) for all n > 1 at finite magnetic field in the 
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vicinity of the given if -point. Again, this contrasts to 
graphene where A and B sublattices carry equal proba- 
bility density for the non-zero LLs. 

At zero energy two contributions add to the density of 
states: First, the zero-energy LL of degeneracy equal to 
the number of flux quanta penetrating the lattice which 
equals the degeneracy of all other LLs. The correspond- 
ing wave function, 

V>o = (0,0, M T 

resembles the wave function of graphene's zeroth LL. On 
the other hand, there is the "non-magnetic" pseudospin 
S p = topological level whose energy remains unaffected 
when a magnetic field is switched on. Its degeneracy 
equals the number of elementary cells in the lattice and 
the corresponding wave functions read 



top 



y/2n + 1 



/ y/n + l (f> n -l\ 





n > 1 



VI. MAGNETIC BARRIER 

From graphene is known 7 that magnetic barriers can 
confine Dirac-Weyl fermions with zero transparency, con- 
trary to electrostatic barriers. We consider a barrier 
parallel to the y-axis of thickness d given by a space- 
dependent magnetic field perpendicular to the lattice 
plane, 



B z (x,y) = B Q(d/2-\x\) , 



(23) 



described by a corresponding vector potential A{x) — 
{0,A y (x),0}, 

{-d/2 x < -d/2 
x for |a;| < d/2 ■ 
d/2 x > d/2 

Stationary states of the Hamiltonian (1) of the form (3) 
satisfy the set of coupled differential equations 



V2£tfj A 
V2£ip H 

V2£?Pb 



id x ijj n -i[ky+ —^A v (x) Vh, (24a) 



-id x ip A + i yk y + —^A y {x)j ipA 

( e£ 2 \ 
■id x ipn + i{k y + -^-A y (x)jip B , (24b) 

/ e f \ 
-id x ip n +ilk y + ^A y (x)j^, (24c) 



where k y is the transverse momentum. Here and in 
the following we have rescaled all quantities by using 



£g = ^c/(eB ) as the unit of length and Vf/£b — u c /y/2 
as the unit of energy, £ = E / {vp / £b)- The wave func- 
tion left of the barrier is of the form (9) with k y — 
£ sin(</?) + d/2, where ip is again the angle of incidence. 
Conservation of momentum parallel to the barrier results 
in the condition 



sin(y) + d/£ = sin(a) 



(25) 



so that the wave function on the right side of the barrier 
reads 



ipm(x) = t 



V2s 

Ve-y 



(26) 



with emergence angle a and final momentum in the IE- 
direction q x = sin(a)£. As in the case of graphene, 7 the 
condition of momentum conservation in the y-direction 
(25) cannot be fulfilled for sufficiently thick barriers, d > 
2£. Finally, the wave function in the barrier region is 
given by 



Uv v+1 (£)-v v (0\ 
V -A, (0 / 
+b ev-vGZ) 



(27) 



where the T> v are parabolic cylinder functions taken at 
argument £ = \/2{k y + x) and r\ = (£ 2 + l)/2. 

The transmission probability T — \t\ 2 can be calcu- 
lated straight-forwardly by applying the boundary con- 
ditions (4). The exact result is lengthy and given in 
Appendix B; it is plotted in Fig. 6. The main features 
of this result can be understood within a semiclassical 
picture, according to which electrons perform cyclotron 
orbits of radius vq = (E/v-p)£ B , inversely proportional 
to the strength of the magnetic field. 26 The semiclassi- 
cal approximation is justified for vq £g. Figure 7(a) 
illustrates three possible cases, (i) For an incident an- 
gle ipi > </? C rit larger than a critical angle the particle 
reaches the opposite boundary and leaves with transmis- 
sion probability T — 1. (ii) At the value of the critical 
angle, ip2 — <p C Tit, the particle leaves the barrier region 
parallel to the barrier, (iii) For angles 953 < </5 C rit, the 
particle makes the full turn under the barrier and is per- 
fectly reflected, T = 0. In this semiclassical picture the 
critical angle ip cr it is obtained from 



dmax = £(1 + sia(v)) 



(28) 

which relates the maximal barrier width <i max that can 
be overcome at incident angle tp and energy £. In fact, 
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FIG. 6: (Color online). The transmission probability T as a 
function of the injection angle for a magnetic barrier of width 
d and at the energy £ for the T3 lattice. In panel (a) the 
parameters are £ = 4.5, d = 1 (red solid line), d = 3 (blue 
dot-dashed line), d = 6 (green dashed line), and d = 8 (yellow 
dotted line). In panel (b) the parameters are d — 2.2, £ — 1.6 
(red-solid line), £ = 2.5 (blue dotted line), and £ = 5 (green 
dot-dashed line). 



(a) 



1.0005 



B=0 B#0 B=0 g 0.9995 h 

E^T 0. 1 
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1 P 2 





FIG. 7: (Color online), (a) Sketch of possible semiclassical 
trajectories of particles entering the barrier region at different 
incident angles, as explained in the text. Panels (b) and (c): 
transmission probability as a function of the normalized width 
of the barrier d/d ma x for the 75 lattice (solid red lines) and 
HCL (dashed blue lines). Panel (b): £ = 6.9 and cp — — 7r/4. 
Panel (c): £ = 3.9 and tp — 0. 



transparent than the HCL, so that the latter is seen to 
develop quantum interference oscillations as a function 
of the barrier width. 



VII. SUMMARY 

We have analyzed transmission properties of barriers 
on the T3 lattice in the low energy approximation using 
the Dirac-like Hamiltonian (1). First we derived and dis- 
cussed the boundary conditions for the wave function at 
a given interface. Wave function components of rim sites 
can be discontinuous, whereas the hub component is con- 
tinuous; this can be well understood in terms of the prob- 
ability currents perpendicular and parallel to the barrier. 
An electrostatic barrier on T3 has enhanced transparency 
compared to the HCL, which is a direct consequence of 
the specific boundary conditions of the 5 = 1 Hamil- 
tonian. When the energy of the incoming wave equals 
half the barrier hight, we even obtain perfect transmis- 
sion independently of the incident angle. This "super 
Klein-tunneling" makes it possible to use a rip-junction 
in order to design a perfect focusing lens without loss; 
this must be contrasted to the 50% of loss of HCL in the 
same configuration. Therefore, we believe that T3 can 
be interesting from the point of view of electron focusing 
and photonic crystals. 27 

Furthermore, we have confirmed numerically the an- 
alytical predictions about the super Klein-tunneling via 
a lattice Green's function method. In addition, we have 
investigated numerically the influence of the intervalley 
scattering in the T3 lattice compared to the HCL and 
identified the crucial role of the flat band. 

For a uniform magnetic field we have discussed the 
LLs; their energies differ from non-relativistic and from 
relativistic systems known so far. T3 exhibits a high 
density of states at zero energy due to contributions from 
a zeroth LL and a topological level. 

Finally, we have investigated the transmission through 
a magnetic barrier in T3 which is found to be similar 
compared to the case of the HCL. Both cases can be 
understood in a simple semiclassical picture from circular 
cyclotron orbits under the barrier. Magnetic barriers of 
sufficient width can be used to confine particles in the 
T3 lattice. 



the curves obtained for the HCL, taking the results of 
Ref. [7] , are almost identical when plotted on the scale of 
Fig. 6, confirming the semiclassical point of view also for 
graphene. 

In Figs. 7(b) and 7(c) we show the transmission prob- 
abilities for both lattices as a function of the normalized 
magnetic barrier width c?/d max at different angles and 
energies where the transmission for d < d max is very 
close to unity. However, the 7~3 lattice reveals as more 
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Appendix A: On the boundary conditions 

We can obtain the boundary conditions from the low 
energy wave equation (1) together with the conserva- 
tion of the probability current. The former states that 
-id t ^ = HV, the latter that d t \ijj\ 2 = -V • j. 

For the case of graphene and with the two component 
spinor ^ — ('0AjV'b) T we get the following condition: 

d t \y\ 2 = (<9 t #)#* + (d t tf*)tf (Ala) 
= -2v F d x Re[r A M ~ 2v F d v Im[r A ^B} (Alb) 

= ~(d X j X + 8yjy) . (A1C) 

Therefore, we can identify the two components of the 
probability current 

j = 2t- F (Re[^ B ],in#XVB]) T (A2) 

At an interface we require the continuity of the prob- 
ability current j x perpendicular to the boundary. This 
implies for the HCL the continuity of the two components 
of the wave function ^ at the interface. 

On the other hand, for the case of the T3 lattice and 
the three-component spinor $ = (vpA, V'h, *Pb) T , we get 

9 t |*| 2 = (<9 t #)#* + («9 t #*)* (A3a) 
= -^F^Re^&^A + V>b)] 

+ V2v F d y Im[^(^ A -4, B )} (A3b) 

= -(d x j x + dyjy) . (A3c) 

In this case the probability current is defined by 
j = V2v F (Re[^ A + ^B)],Im[^(^A - ^b)]) T (A4) 
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